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Abstract 

We construct the creation operator of a vortex using the methods developed for 
monopoles. The vacuum expectation value of this operator is interpreted as a dis- 
order parameter describing vortex condensation and is studied numerically on a 
lattice in SU(2) gauge theory. The result is that vortices behave in the vacuum in 
a similar way to monopoles. The disorder parameter is different from zero in the 
confined phase, and vanishes at the deconfining phase transition. We discuss this 
behaviour in terms of symmetry. Correlation functions of the vortex creation oper- 
ator at zero temperature are also investigated. A comparison is made with related 
results by other groups. 
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1 Introduction 



Recently there has been renewed interest in understanding the role of vortices 
in the mechanism of colour confinement [1-3]. Vortices were originally intro- 
duced in the continuum theory as string-like topological defects [4], and have 
been studied both in the continuum and on the lattice [5]. In the notation of 
Ref. [4], a vortex creation operator B(C), for a gauge group SU(N), can be 
associated to each closed oriented curve C and has the following commutation 
relation with a Wilson loop W(C), bounded by a closed curve C: 

W{C) B{C) = B(C) W(C) exp(2irin/N) (1) 
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where n is the winding number of C around C. In SU(2) such a vortex 
contributes a factor (—1) to each Wilson loop with an odd linking number 
with the vortex line. On a lattice, the world-sheet M of the vortex can be 
associated to a surface on the dual lattice. A vortex creation operator is defined 
by twisting the plaquettes in the action, whose duals belong to that surface. 

While creating a vortex is a well defined procedure, detecting vortices in lat- 
tice gauge configurations is more difficult. Attempts in this direction stemmed 
from [6], and results were obtained by several groups [7], indicating that vor- 
tices can play a role in confinement. However, it is not clear that this procedure 
detects the vortices as defined in Eq. 1 (see e.g. [8,9]). The current picture of 
the role of vortices in colour confinement is based on Eq. 1: at low temperature, 
vortices disorder the Wilson loop to produce the well-known area law; at high 
temperature, vortices are suppressed and the Wilson loop obeys a perimeter 
law. 

In three dimensions, a conserved topological quantum number can be associ- 
ated to a vortex, and a vortex creation operator is defined as a local complex 
scalar field. The dynamics of this scalar field is described by an effective La- 
grangian with a global dual symmetry [4] . The breaking of this symmetry 
is responsible for the different phases of the theory. The scalar field 4>(xq) in 
three dimensions has a non-trivial commutation relations, analogous to those 
of Eq. 1, with any Wilson loop encircling the point x . Operators that cre- 
ate such topological excitation are known in statistical mechanics as "disorder 
operators" [10] and the expectation value of a disorder operator is a disor- 
der parameter and is related to the free energy of the associated topological 
excitation. 

In four dimensions the situation is less clear. The vortex corresponds to a 
string-like topological defect and the dual symmetry does not emerge natu- 
rally. In this respect, there is a fundamental difference between the vortex 
and the monopole picture of colour confinement in terms of symmetry: for the 
monopoles, a topologically conserved current exists, whose zero component 
is the generator of the dual symmetry, and a picture based on a mechanism, 
namely the dual Meissner effect [11], emerges as a consequence of the breaking 
of the dual symmetry [12]. For the vortices the symmetry pattern is not under- 
stood. In close analogy with previous work on abelian-projected monopoles, it 
is nonetheless possible to define a vortex creation operator obeying the above 
commutation relation, which can be associated to the free energy of the vor- 
tex. Starting from this point of view, in this paper we study the role of vortices 
across the deconfinement phase transition, using the techniques developed in 
Refs. [13,14,12] for investigating the condensation of monopoles defined via 
the abelian projection. 
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A vortex creation operator /i is introduced as a disorder operator for the 
SU(2) lattice gauge theory. When applied to a gauge configuration, /i creates 
a vortex associated to a string closing in the z direction via periodic boundary 
conditions (PBC). As well as for other systems, both in statistical mechan- 
ics [10,15,16] and in quantum field theory [17,13,14,12], such an operator can 
be constructed without any reference to the dual theory. In Sect. 2 we give the 
explicit construction of this operator. While our study was in progress, other 
papers appeared which addressed similar problems by using quantities related 
to the vortex creation operator examined in this work. The explicit relation- 
ship between the operator studied here and other quantities in the literature 
is reviewed in Sect. 3. 

Correlation functions of the vortex creation operator at zero and finite temper- 
ature yield useful informations on the behaviour of these topological defects 
in the SU (N) Yang-Mills vacuum. As in previous studies on monopoles, it is 
convenient to compute 



and analogous quantities for n-point correlators. Unlike /i, p can be determined 
with good accuracy and provides all the needed information. The vacuum ex- 
pectation value (vev) of /i proves to be a disorder parameter for the finite 
temperature deconfining phase transition, as reported in Sect. 4. Another in- 
teresting quantity to be investigated is the correlator (/i(t )/i(t + t)) at zero 
temperature. 

This approach can be generalised to SU (N) gauge group and to vortex loops 
of generic geometry, also not wrapping through PBC. Work is in progress on 
these aspects. 



2 Vortex creation operator 

The vortex creation operator \x is constructed by the same technique used 
for the monopole creation operator in previous works [13,14]. The vacuum 
expectation value of \i is defined as the ratio of two partition functions: 



d 



log(/x) 



(2) 






where Z is the usual partition function with the Wilson action, 




(4) 
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and S is obtained from S by twisting a line of plaquettes in the Oy plane: 

P 0y (t , x > x , y , z) h-> -P 0y (t , x > x , y , z), Wz (5) 

The net effect of this transformation is best understood by performing a series 
of changes of variables in the functional integral as follows. The first change: 

U y (t + l,x > x ,y ,z) h-> -U y (t + l,x > x ,y ,z), Wz (6) 

yields 

(^t ,x ,y )) = Z- 1 J[dU]e-^ W M ( 7 ) 

where is defined by the following transformation of the Wilson action: 

Poy(t = t + l,x> x , y , z) ^ -P 0y (t = t + l,x> x , y , z), Wz 

Pox ^ Pox: Pqz i— > Pqz 
* Pxy(t = t + 1, x , y , z) ^ -P xy {t = t + 1, x Q , i/o, z), Wz 
P xy (t = t + l,N a -l, y , z) ^ -P xy (t = t + l,N a - 1, yo, z), Wz 




z 



Fig. 1. Set of plaquettes changing sign in the hyper-plane io + 1 after the first change 
of variables described in Eq. 6. The plaquettes in the plane xy correspond to the 
location of the vortex lines at (xo,yo), and at (N s — l,yo) due to PBC in the x 
direction. 

The change of sign in the P 0y plaquettes introduced at t — t in Eq. 5 has 
been shifted at t — to + 1 and two vortices have been created in the (x, y) 
plane for t = t + 1 at (x , y ) and (N s — l,yo). It is worthwhile to remark that 
the second vortex at (N s — l,yo) is due to the fact that we are working in a 
finite volume with PBC in the x direction. In an infinite volume, or with free 
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BC, the vortex at (xo,yo) would have been created alone. In the z direction, 
the vortices either extend to infinity (for an infinite lattice), or they form a 
closed loop due to PBC. The set of xy plaquettes changing sign after this first 
change of variables is shown in Fig. I. 

The same change of variables can be iterated at successive times, yielding, 
after n iterations: 

(»(to,x ,y )) = Z' 1 J[dU]e-^ (n) M ( 8 ) 
and is obtained from the Wilson action via: 

Po y (t = t + n,x > x , y , z) h-> -Po y (t = t + n,x> x , y , z), \/z 

P xy (t = t + l,x , y , z) ^ -P xy (t = t + l, x , y , z), \/z 

P xy (t = t + l,N 8 -l, y , z) ^ -P xy (t = t + l,N 8 - 1, y , z), \/z 

P xy (t = t + n, x , y , z) i-> -P xy (t = to + n, x , y , z), \/z 
k Pxy(t = t + n,N s - 1, y , z) ^ -P xy (t = t + n,N s - 1, y , z), \/z 

Such a configuration corresponds to the propagation of the two vortices, from 
time t — to to time t = t + n. 

Correlations of \i operators at different times {t±, . . . ,t n } are defined by re- 
peating the transformation in Eq. 5 for each For a two-point correlation: 

T(t) = (/x(t , x , y ) //(t + t, x Q , y )> (9) 

the change in the P 0y plaquettes is reabsorbed when n — t and the correlator 
describes the propagation of a pair of vortices from the time t to the time 
t + t. 

We remark here that these prescriptions create vortex lines that are closed in 
space, due to PBC in the z-direction, and propagate in time. 

At large values of t, T decreases exponentially to an asymptotic value, which, 
by cluster property, is (/i) 2 , the square of the disorder parameter: 

r(t) ~ Ae- mt + (/i) 2 (10) 

At T = 0, (//) can be extracted from the large-t value of the correlator, 
according to Eq. 10. 

At finite temperature, there is no propagation in time, and a direct measure- 
ment of (/t) is needed by use of a single operator. In order to have a consistent 
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implementation of the changes of variables described above, C* boundary con- 
ditions are needed in the time direction, as explained in [14]. 

Closed vortex lines propagating in time can be created in a field configuration 
with the same prescription, but with a different modification of the action. 
For instance, using for the modified action: 

Po y (to, %o < x < xi, y , z < z < zt) i-> -P 0y (t , x < x < xi, y , z < z < z x ) 

(11) 

gives a vortex line encircling the rectangle in the xz plane: 

R= {(x,y,z) : x < x < x 1 ,y = y ,z < z < zi} (12) 
The same definition is used in a recent publication [21]. 

3 Comparison with related works 



Our disorder parameter is strictly related to observables introduced by other 
authors [1-3]. The detailed comparison is as follows. 

The paper [2] presents a computation of the free energy of a vortex pair. The 
latter is defined as: 

where 



Z((3,(3')= [dU]e 



S{P,P) = \[P £ TrP + Z?' £ TrP 
- '•■ p^m PeM 



and T = (A r t a)~ 1 is the temperature of the system. M indicates a set of pla- 
quettes with modified coupling f3' = —j3, which can be seen as the plaquettes 
transversed by the vortex string. The following two cases are examined: 

• a vortex solution is placed at (xq, yo, Zo, 0) and an anti- vortex at (xo, yo, zq + 
d, 0), with a straight string of twisted plaquettes connecting them; 

• a single vortex is placed in the middle of the lattice and the string extends 
to the boundary in the z-direction, where free boundary conditions are used. 

In both cases the modification of the action is done in all time-slices. 
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Using the notation introduced in the previous section: 




(14) 



and now S' 



is obtained from S by: 



Pxy(t, Xq, yo, z 



20 + 1) 



P xy (t,x ,y ,z = z + 1), Vi 



< 



Pxyit, Xq, l/o, z 



z + d) 



Pxy{t, Xq, l/o, Z = Z + d) , W 



With a t ^ z relabelling of the axes, the correlator T defined in this paper by 
Eq. 9 yields the free energy of two vortex pairs as defined in [2] at distance 



At T = we use symmetric lattices (N s = N t ) and measure the time correlator 
of two fi operators at x — y — N s /2, which is in fact the correlator of 
two pairs, when the effect of periodic boundary conditions are taken into 
account. At non-zero temperature, we use a single vortex operator, which 
actually introduces two vortices, again due to PBC. The extra vortices that 
are created by PBC are N s /2 lattice spacings away from their partners. If 
N s ^> £, where £ is the correlation length, then our disorder parameter (/i) is 
the square of exp (-F/T). 

The authors of Ref. [3] are concerned with the behaviour of the 't Hooft 
loop [4] in hot QCD, and relate the dual string tension to the wall tension 
of Z N domain walls. For SU(2), the lattice definition of the loop operator is 
given by [18]: 



where S is the surface in the xy plane bounded by the curve C. The vev of the 
operator defined in Eq. 15 is again the ratio of two partition functions: one 
with a modified action divided by the (standard) Wilson action. The modified 
action is obtained by changing the sign of the plaquettes in S. The definition 
given in sect. 2 coincides exactly with Eq. 15, for the curve C depicted in Fig. 2, 
after a relabelling of the axes y ^ z. It is easy to realise in this formulation that 
our vortex creation operator is precisely a 't Hooft loop operator: every Wilson 
loop with non-trivial linking to the curve C receives a (—1) contribution from 
the vortex. 

Finally, the authors of Ref. [1] compute the free energy of a closed 577(2) vor- 



N s /2. 




(15) 
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z 




Fig. 2. The operator p, introduced in sect. 2 corresponds to V[C) for the curve C 
depicted here. The box represents a time-slice of the whole lattice. 

tex, which is defined from the logarithm of the ratio of two partition functions: 

exp{-F(r M ,)} = ^ (16) 

Zij^) is defined as usual by multiplying a given (co-closed) set V^ u of plaque- 
ttes in the action with an element t m „ of the center of the gauge group. For 
SU(2), the only non-trivial possibility is — — 1. For the set of plaquettes: 

V»v = {P»v{t, x, y, z) : (n, v) = (0, 2), t = t Q , x < x < x u y = y , z < z < z x } 

the definition in [1] coincides with the one in Eq. 11. 

4 Numerical results 

In this paper, we present data only for open vortex lines wrapping in the z 
direction by PBC. Work is in progress to study closed vortices as done in [21]. 

Due to the exponential in its definition, \x has large fluctuations, which make a 
direct measurement a challenging task. Instead, as in previous studies [13,12], 
we define 

P = Jp lo s(^ (*o, %o, yo)) = (S)s ~ (S)s (17) 

Being the difference of the expectation values of two actions, p can be easily 
computed, and proves to yield all the relevant information. At finite tem- 
perature, p is expected to have a sharp negative peak in the critical re- 
gion [13,14,12], if (p,(to,xo,yo)) is related to the deconfinement phase tran- 
sition. 
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Fig. 3. Comparison of p defined with the vortex and the monopole creation operator. 

Our results for p are displayed in Fig. 3, for several lattice sizes. For com- 
parison, we also report a plot of the corresponding quantity for the vev of a 
monopole creation operator on a 12 4 x 4 lattice [12]. Already at a qualitative 
level, it is clear that, at weak coupling, p is negative, with its absolute value 
increasing with the volume. This behaviour implies that (p) vanishes for large 
(3 in the thermodynamical limit. At low (3, p is compatible with for all the 
volumes considered in this study, which means that (p) has a non-zero value 
in the infinite volume limit in the confined phase. The presence of the nega- 
tive peak connecting the confined and deconfmed phases suggests that vortices 
play a role at the deconfinement transition. This behaviour is in agreement 
with the results in [1,2]. As in the case of monopoles [12], we can perform a 
finite size scaling analysis of the data presented above. In neighborhood of the 
critical coupling /3 C , if (p) is a disorder parameter for the deconfming phase 
transition, in the infinite volume limit we have 



(p)^(p c -py 



being 5 the corresponding critical index. In a finite volume, the previous equa- 
tion is replaced by 



(p) = (P c -P) 5 $(Ay£) (19) 
where £ is the correlation length and $ is a function of the ratio N s /£. Since 



(20) 
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Eq. (19) can be written as 

(ri = L-""*{L 1 '"(p e -0)) (21) 
being v the critical index of the correlation length. 

Eq. (21) implies 

jfc = f (L 1/v (& - (3)) (22) 

i.e. the ratio p/L 1 ^ is an universal function of the scaling variable 

x = L 1 '" (& - p) 

o i . , . 1 . 1 . 1 



-0.5 



J -1 - 



-1.5 




By guessing that [12] 
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Fig. 4. Plot of rescaled data. 



15 



5 

— + c 

X 



(23) 



it is possible to extract from our data the critical exponents S and v and the 
critical coupling (3 C . We perform three different fits: a fit with three parameters, 
a fit to 5 and v at fixed (3 C , using the value in [19], and again a two parameter 
fit to 5 and v at fixed /3 C , using the value obtained in the three parameter 
fit. The error is estimated from the variation in the fitted values when using 
the different methods described above and when the points corresponding to 
smaller correlation lengths are excluded. Our results are: 
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/3 C = 2.30(1) 
5 = 0.5(1) 
v = 0.7(1) 

The values of (3 C and z/ are in good agreement with the values obtained in 
[12], while the value of 5 varies by two standard deviations, when compared 
to previous results. 

Fig. 4 shows how well the scaling relation is obeyed with the values of (3 C and 
v from the above fit. 

Assuming that in the weak coupling limit all link variables are close to unity, 
and remembering that the twisted action is obtained by flipping the sign of 
2 x N t x L plaquettes, a naive prediction would yield: 

p~-16L + const (24) 

This is a justification for the linear dependence used in the ansatz above. 
However, the actual value of the coefficient a does not need to be the one in 
Eq. 24, since the configuration with all links set to the identity is not the true 
minimum of the modified action which enters the definition of p. 

To check this behaviour, we have measured p for different volumes and very 
large values of the coupling. In this region, we expect the data to be indepen- 
dent of /3. Within the errors the data do lie on a straight line as predicted by 
Eq. 24. The fact that the slope is negative ensures that the disorder operator 
vanishes in the weak coupling phase in the thermodynamic limit. 
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The behaviour of p at low /3's is displayed in Fig. 6: the value at low (5 remains 
bounded from below and consistent with zero for increasing volumes, which 
guarantees that p does not vanish in the thermodynamic limit above the phase 
transition. 
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Fig. 6. p vs. for different spatial volumes. 

Another quantity that can give information on the behaviour of center vortices 
in the vacuum of SU(2) is the correlation function of two p operators as 
a function of distance at zero temperature. When relabelling the axes, as 
discussed in the previous Section, this quantity is identical to the free energy 
recently discussed in [2]. 

Again, our numerical computation is performed in terms of p2, now defined as 



P2{t) = \og(fi (t , x , y ) p (t + t, x Q , y )) 



(25) 



which is the difference between two actions. Since the usual Wilson contribu- 
tion is ^-independent, the whole dependence on t is brought by S. 

In Fig. 7 we display the behaviour of S on a 16 4 lattice at j3 = 2.5, where the 
system is confined and scaling is supposed to work. 

The t dependence of the correlation is fitted using two different fitting func- 
tions, both depending on three parameters (a, b, m) : 



-rat 



p 2 (t)=a- 
p 2 (t) = ae 



+ b 



t 

-mt + b 
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Fig. 7. The modified action S per unit of voiume as a function of the distance 
between the vortex and the antivortex on a lattice 16 4 at j3 = 2.5. 

these two fits will be named respectively fit I and fit II in what follows. The 
functional dependence in fit I describes a Yukawa potential between vortices, as 
studied in Ref. [2]. In both cases, b is the asymptotic value of p 2 and is expected 
to be different from zero if vortices are condensed. Both ansatz fit the data 
relatively well, despite the fact that the masses obtained are quite different. 
The values of x 2 & ve an indication that fit II could be a better description 
of the data, although it is impossible to draw any robust conclusion given 
that it is very difficult to disentangle the logarithmic correction in t which 
differentiate fit I from fit II. The outcome of the fits is summarised in Table 1 
and the two fitting curves are superimposed on the data in Fig. 7. The value 
of the constant b is determined quite accurately due to the plateau in the data 
points at large t. It is interesting to remark that a turns out to be negative 
in both cases up to 90% CL, while the error on the fitted mass is very large. 
This is mainly due to the fact that large variations in m can be reabsorbed 
by variations in a in the range where data points are available. 



fit type 


a 


b 


m 


X 2 /d.o.f. 


fit I 


-0.00018 


0.34923 


0.650 


2.6 


fit II 


-0.00035 


0.34923 


1.302 


1.8 



Table 1 

Fit results for the correlation of vortex operators 

We notice from fit II that the correlation length is around 1.3 lattice spacings. 
This justifies the view of the pair of vortices (the one at the center of the 
lattice and the one at the border produced by PBC) as independent. 
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5 Discussion and future outlooks 



In previous works [12], we produced evidence of dual superconductivity of the 
QCD vacuum, supporting the mechanism of Refs. [11]. A peculiar feature of 
this phenomenon is that monopoles defined by different abelian projections 
all condense [12]. This possibility was suggested in Ref. [22] and lattice data 
support it, showing that confinement is related to condensation of magnetic 
charges defined by a few abelian projections. 

If all or a large class of abelian projections are equivalent, there are infinitely 
many physically equivalent disorder symmetries. We already observed that 
this fact is not inconsistent, but suggests that maybe a more fundamental 
dual symmetry pattern exists, which can manifest itself as condensation of all 
these magnetic charges [23]. 

Now it is found that also vortices show condensation. We were not able to 
associate a dual conserved topological quantity to vortices in 3+1 dimensions, 
contrary to what happens in 2+1 dimensions. However we consider what is 
observed here an important information on the way to understand the features 
of the dual description. 

We are extending the analysis to the vortices of the SU(3) gauge theory 
and to some aspects of closed vortex lines. We are also trying to understand 
what happens in the presence of quarks. Preliminary data on condensation 
of monopoles, as defined by the abelian projections, show that the presence 
of quarks does not affect the dual superconductor picture. This is in agree- 
ment with the idea that the gross features of the QCD vacuum, including the 
mechanism of confinement, are determined already at N c = oo [24]. 

We are now trying to understand how to extend to full QCD the analysis of 
vortices. 
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